Electron Transport Through Molecules: Gate Induced Polarization and Potential Shift 
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We analyze the effect of a gate on the conductance of molecules by separately evaluating the gate- 
induced polarization and the potential shift of the molecule relative to the leads. The calculations use 
ab initio density functional theory combined with a Green function method for electron transport. 
For a general view, we study several systems: (1) atomic chains of C or Al sandwiched between Al 
electrodes, (2) a benzene molecule between An leads, and (3) (9,0) and (5,5) carbon nanotubes. We 
find that the polarization effect is small because of screening, while the effect of the potential shift 
is significant, providing a mechanism for single-molecule transistors. 

PACS numbers: 73.40.Cg, 72.10.-d, 85.65.-|-h 
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Electron transport through molecules sandwiched be- 
tween two metallic electrodes has been attracting increas- 
ing attention both for fundamental reasons and because 
it may form the basis of a future molecular electronics 
technology.^-^ To control the transport and to realize 
single-molecule-based transistors, a gate is usually ap- 
plied to the two-terminal lead- molecule- lead (LML) sys- 
tem. This gate terminal can be (1) global - i.e. much 
larger than the molecular region - made by constructing a 
plate capacitor near the devicef^*^*^ or (2) local, grown, 
for instance, by using chemical vapor deposition beneath 
the device regioni or simply realized by the sharp tip of 
a scanning probe microscope i^ 

The application of a gate to a LML system has two 
conceptually distinct consequences. First, the voltage on 
the gate will create a strong electric field in the direction 
perpendicular to the transport. Second, because of the 
field component along the transport direction, there will 
be a potential shift in the molecular region with respect 
to the leads. Both effects may influence significantly the 
electron transport through the molecule. 

Theoretically, a tiny local gate may be modeled in ab 
initio calculation by a finite flat surface with a constant 
electrostatic potential. Previous ab initio calculations 
have evaluated the polarization effect of a local gate for 
one caseS. by puting the molecule into a plate capacitor, 
and the total effect of a local gate in two othersiSiii by 
solving the Poisson equation for the molecular system in- 
cluding the gate. However, in most recent experiments a 
golbal gate is used, and its polarization electric field can 
be quite different from that of a tiny local gate. Addi- 
tionally, the issue of which effect of the gate is dominant 
has not been addressed, either experimentally or theo- 
retically. Thus, the individual contributions of the two 
effects remain unclear. 

In this paper, we investigate these two effects - po- 
larization and potential shift of a global or local gate - 
separately from first principles. The approach is through 
the well-known combination of density functional the- 
ory (DFT)^ for the electronic structure with the non- 
equilibrium Green function (NEGF) methodi2*i^ for elec- 
tron transport. Our implementation was developed pre- 



viously for two terminal systems^ We use periodic 
boundary conditions (PBC) for the DFT calculation, as 
a result, the geometry of the LML system is accurate 
without any artificial surface effects. Using this method, 
we calculate the molecular conductance under an exter- 
nal polarization electric field or a potential shift in the 
molecular region, both induced by a gate. 

To arrive at a general view of gate effects, here we 
investigate three kinds of typical systems: (I) atomic 
chains, C and Al, sandwiched between two Al(OOl)- 
2V2 X 2\/2 leads; (2) a benzene molecule between two 
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FIG. 1: Structures of the LML systems calculated: (a) car- 
bon atomic chain sandwiched between Al(OOl) leads, (b) alu- 
minium chain between Al(OOl) leads, (c) benzene molecule 
connected to Au(OOl) leads through S atoms, (d) (9,0) car- 
bon nanotube, and (e) (5,5) carbon nanotube. The dotted line 
indicates the interface between the device region (extended 
molecule) and the leads. The gray bar labeled Vg marks the 
region where the potential is shifted by the gate. 



Au(001)-2-\/2 X 2-\/2 leads; (3) carbon nanotubes, both 
semiconducting (9,0) and metallic (5,5). The atomic 
structures of these LML systems are shown in Fig. 1. 
Vg denotes the local gate potential shift or the locally 
projected potential shift induced by a global gate. 

For zero bias, the density matrix D of the device region 
is simply related to the Green function by 
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where Gjj (E) is the retarded Green function of the de- 
vice region, / is the Fermi function, and /i^ the chem- 
ical potential of the leads. On the other hand, Gd{E) 
is determined by the self-energies for the two leads to- 
gether with the molecular Hamiltonian (H^i) given by 
DFT based on the density implied by D^i. After these 
two steps converge selfconsistently, the transmission co- 
efficient at any energy, T{E,Vb), is calculated from the 
Green function. The conductance, G, then follows from 
a Landauer-type relation. 

We use a uniform electric field E to model the field per- 
pendicular to the transport direction caused by a global 
gate. Because we use periodic boundary conditions for 
the DFT calculatioufi^ there will be a potential jump 
at the boundary; this has, however, no unphysical effect 
since it is in the deep vacuum. Under the field there is 
an additional term, — r • E, in Hd- In the self-energies 
for the semi- infinite leads, the same uniform electric field 
is applied. As a result, the whole infinite LML system 
will be under a uniform field. Then, the self-consistency 
is carried out just as in the zero field case. 

To simulate the potential shift along the transport di- 
rection induced by a gate, we apply a potential shift di- 
rectly to the molecular region (M), 
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where S^i is the overlap matrix of the device region. Be- 
cause the shift is applied to a matrix element when ei- 
ther orbital index is in M, the potential shift will be 
slightly smeared around the molecule-lead contact; we 
believe this is more realistic than a step-function change. 
This potential shift may be applied in a non-selfconsistent 
(non-SC) or selfconsistent (SC) way. In the non-SC 
scheme the shift is added to the initially converged Hjj 
once without further iterations. In the SC scheme AH/^ 
is added to the initially converged H/^ for each iteration 
until Hd and T)d is converged. 

A gate potential shift will induce charge accumula- 
tion in the molecular region, AQmi given by AQm = 
Timi^D • Djj]. Here Tr^ means trace over only the 
molecule part, and D/) is determined by the initially con- 
verged Hd plus AH^i according to Eq. QJ . This charge 
accumulation in the molecular region will induce opposite 
charge accumulation on the lead surfaces (AQl)- 

The difference between the two schemes is that AQl 
is not taken into account in the non-SC scheme while it 
is included in the SC scheme by the charge normalization 
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FIG. 2; Equilibrium conductance as a function of a uniform 
electric field applied to (a) the C7 and AI7 atomic chains, and 
(b) the benzene molecule anchored by S atoms. The direction 
of the field with regard to the plane of the benzene ring is 
indicated in (b). Note the very small variation. 



process in each iteration. The SC scheme, therefore, pre- 
serves the charge neutrality of the device region. Note 
that in the non-SC scheme AQm is not compensated at 
all while in the SC scheme AQm is compensated fully by 
AQl- (For the SC scheme, we have to incorporate larger 
parts of the leads into the device region so that AQl can 
be accommodated.) The real situation obviously lies be- 
tween these two extremes: AQm is compensated partly 
by AQl and partly by the charge accumulated on the 
gate. To reveal the difference between these two limits 
and to check that our approach is reasonable, we carry 
out calculations with both schemes. 

We adopt a LCAO-like numerical basis set to ex- 
pand the wave functions^ and make use of op- 
timized TrouUier-Martins pseudopotentialsi^^ for the 
atomic cores. For the benzene system with gold leads, a 
high-level double zeta plus polarization basis set is used 
for all atoms, while for the systems without gold a single 
zeta basis set is used. The PBE version of the general- 
ized gradient approximation^" is adopted for the electron 
exchange and correlation. 

In Fig. 2 we show the calculated equilibrium conduc- 
tance as a function of the global electric field for the 
atomic chain and benzene systems. Note that the field 
applied is very large, approaching the breakdown limit. 
Among the atomic chains, the polarization is larger for 
C than for Al because the screening is stronger in the 
latter. For the benzene molecule, the polarization effect 
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FIG. 3: Transmission functions of the (9,0) and (5,5) nan- 
otubes under zero and 0.5 V/A electric fields, 'uniform field' 
('local field') means that the whole nanotube (a finite part of 
the nanotube) is under the field. 'coupling_0' and 'coupling_l' 
indicate the different ways to treat the coupling bwteen the 
regions with and without the field, as explained in the text. 
Note that the effect is very small around the Fermi energy. 



is stronger when the field is perpendicular to the ring 
than when it is parallel, which is related to the charac- 
teristics of the large tt bond. However, as can be seen, 
the overall polarization effect is negligible in all cases, 
indicating that screening plays a very important role in 
electron transport under a perpendicular electric field. 

In Fig. 3 we show the equilibrium transmission func- 
tions of the (9,0) and (5,5) carbon nanotubes under zero 
and 0.5 V/A electric field. The field is applied either 
uniformly (as described above) or locally. A local field 
means that it is applied to only part of the device (in- 
dicated by the gray bar in Fig. 1) and is not included 
in the lead self-energies. We treat the coupling between 
the regions with and without the field in two approximate 
ways: with or without the field, denoted 'coupling_0' and 
' coupling_l ' , respectively. 

For zero field, T{E) shows clear steps because of the 
quantized electronic states in the directions perpendic- 
ular to the transport. For the metallic (5,5) tube, the 
equilibrium conductance is 2 conductance quanta (Go ~ 
2e^//i) - as is well known, two subbands cross the Fermi 
energy in this case. For the (9,0) tube, there is a small 
gap in T{E) around the Fermi energy - it becomes semi- 



conducting because of the curvature of the tube. Under 
the very strong field of 0.5 V/A, T{E) still has step struc- 
ture and is, in fact, almost the same around the Fermi 
energy. If this strong field is applied locally, however, 
then the step structure is destroyed, but T{E) is still al- 
most the same around the Fermi energy. Thus, the effect 
of polarization is only felt far from the Fermi energy, in- 
dicating that screening plays a very important role just 
as in the atomic-chain and benzene systems. 

From the similarity of our results for completely differ- 
ent systems we come to a general conclusion: screening is 
significant for electron transport through nano- junctions, 
and, as a result, the effect of a polarizing electric field on 
transport is small. 

Fig. 4 shows the equilibrium conductance G as a func- 
tion of gate potential shift for four LML systems, cal- 
culated by the SC and non-SC schemes. Also shown is 
the charge accumulation in the molecular region, AQm, 
due to the gate potential shift. The first thing to no- 
tice is that the equilibrium conductances of all four sys- 
tems are significantly modified by the gate potential shift. 
For some systems, like the Al chain or benzene molecule, 
G can vary sharply from near zero to more than one 
due to only a small change in the gate potential. This 
behavior provides a good mechanism for single-molecule 
transistors. For benzene, the variation is very sharp just 
around Vg = 0; this is related to a resonance near the 
Fermi energy caused by the two additional Au atoms at 
the contactsSi [see the contact structure in Fig. 1 (c)]. 

Different potential shifts lead, of course, to different 
AQm- The function AQm{E) shows a sort of step-like 
structure (Fig. 4). In contrast to the Coulomb blockade 
in which the net charge shows sharp steps, here the steps 
are smeared out because of the broadening of the molec- 
ular orbitals due to the strong molecule-lead coupling. 

As mentioned previously, the difference between the 
non-SC and SC approaches is that AQl is not included 
in the former while probably over included in the latter. 
However, it turns out (Fig. 4) that the SC and non-SC 
results are actually in good agreement for all four sys- 
tems if the gate potential shift is not too large, | Vg | < 
2 V. Note that for a global gate the potential shift here 
corresponds to the local projection of the global potential 
shift, which means that a small | Vg \ can correspond to 
a much larger global gate potential. Even for these rela- 
tively small potential shifts, AQm is already significant, 
implying that AQl should be as well. The good agree- 
ment between the SC and non-SC results indicates that 
the effect of this charge accumulation and the resulting 
electric field is unimportant for the electron transport. 
This conclusion is similar to that for the external elec- 
tric field above. This result also justifies that the present 
SC and non-SC results are reasonable because the real 
situation lies somewhere between these two extremes. 

In summary, we have investigated separately the effects 
of an external electric field and a potential shift induced 
by a gate on the molecular conductance of five different 
lead-molecule-lead systems. We find that the polariza- 
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FIG. 4: Equilibrium conductance as a function of gate potential shift {Vg, a positive value means an upward shift of the electron 
energy) calculated by the non-SC approach (solid line) and SC approach (solid dots), for the four LML systems as indicated. 
The charge accumulation in the molecular region {AQm, in units of electrons) is shown by a dotted line. 



tion effect of the external field is small because of screen- 
ing while the potential shift is significant. The latter 
provides a good mechanism for single-molecule transis- 
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